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Abstract 



We investigate energy-efficiency issues and resource allocation policies for time division multi-access 
(TDMA) over fading channels in the power-limited regime. Supposing that the channels are frequency-flat 
block-fading and transmitters have full or quantized channel state information (CSI), we first minimize 
power under a weighted sum-rate constraint and show that the optimal rate and time allocation policies can 
be obtained by water-filling over realizations of convex envelopes of the minima for cost-reward functions. 
We then address a related minimization under individual rate constraints and derive the optimal allocation 
policies via greedy water-filling. Using water-filling across frequencies and fading states, we also extend 
our results to frequency- selective channels. Our approaches not only provide fundamental power limits 
when each user can support an infinite number of capacity-achieving codebooks, but also yield guidelines 
for practical designs where users can only support a finite number of adaptive modulation and coding 
(AMC) modes with prescribed symbol error probabilities, and also for systems where only discrete-time 
allocations are allowed. 

Keywords: Convex optimization, water-filUng, time division multi-access, fading channel. 

I. Introduction 

With battery operated communicating nodes, energy efficiency has emerged as a critical issue 
in both commercial and tactical radios designed to extend battery lifetime, especially for wireless 
networks of sensors equipped with non-rechargeable batteries. Because the energy required to 
transmit a certain amount of information is an increasing and strictly convex function of the 
transmission rate [1], energy-efficient resource allocation has attracted growing attention [2]-[8]. 
Among them, [2], [3], [4], [5] dealt with energy-efficient designs based on packet arrival and delay 
constraints over additive white Gaussian noise (AWGN) channels; while [6] and [7] considered 
energy-efficient scheduling for time division multi-access (TDMA) networks over fading channels, 
where the data of each user must be transmitted by a given deadline. Recently, [8] minimized 
transmit power of orthogonal frequency-division multiplexing (OFDM) systems using quantized 
channel state information (CSI) of the underlying fading channel. 

Resource allocation for fading channels also remains a popular topic in information theoretic 
studies. However, optimization has been typically carried out to maximize rate (achieve capacity) 
subject to average power constraints. Assuming that both transmitters and receivers have available 
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perfect CSI, Tse and Hanly derived the ergodic capacity [9] as well as the delay-limited capacity 
regions [10] along with the optimal power allocation for fading multi-access channels, while Li 
and Goldsmith found the ergodic [11] and outage capacity regions [12] as well as optimal resource 
allocation policies for code division (CD), time division (TD) and frequency division (FD) fading 
broadcast channels. As regarding delay-limited capacity (a.k.a. "zero-outage capacity"), [13] and 
[14] extended the results of [10], to characterize the outage capacity regions for single-user fading 
channels and multi-access fading channels, respectively. 

In this paper, we re-consider these information theoretic results pertaining to rate efficiency and 
investigate optimal resource allocation for fading channels from an energy efficiency perspective. 
Specifically, we seek to minimize energy/power cost under average rate constraints for TDMA 
fading channels, given perfect or quantized CSI, at the transmit- and receive-ends. As stated in [11], 
TD and FD are equivalent in the sense that they exhibit identical ergodic capacity regions and 
corresponding optimal resource allocation policies. Thus, our results apply also to FDMA fading 
channels. Unlike [2] -[7], we do not impose delay constraints in our energy minimization problems. 

We first study the problem of minimizing total power given a weighted average sum-rate 
constraint for the block flat-fading TDMA channel (Section III). This is dual to [1 1], where rate was 
maximized under a sum-power constraint. Note that we impose a weighted sum-rate constraint for 
the multi-access channels while [11] consider the sum-power constraint for the broadcast channel. 
We then optimize energy-efficiency when each user can only support a finite number of adaptive 
modulation and coding (AMC) modes. The second problem we consider is power minimization 
under individual rate constraints, which is the general case for multi-access channels (Section IV). 
Rate maximization under individual power constraints has been addressed via superposition coding 
and successive decoding in [9] and [10]. Here we formulate and solve energy minimization under 
individual rate constraints for TDMA when users have infinite-codebooks, or, a finite number of 
AMC-modes which requires only quantized CSI to be fed back from the receiver to the transmitters. 
Section VI provides some numerical results, followed by the conclusions of this paper. 
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II. Modeling Preliminaries 



We consider a set of K users linked wirelessly to a single access point and adopt a discrete-time 
multi-access Gaussian channel model as in [9], [14]: 



where Xh{n) and /ife(n) are the transmitted signal and the fading process of the kih user, respectively, 
and z{n) denotes AWGN with variance cr^. Different from [9] and [14], we confine ourselves to 
TDM A where each user transmits in a dedicated time fraction, not overlapping with other users; 
i.e., when Xk{n) 7^ in (1), we have Xi{n) — for \/i ^ k. We also assume that the fading 
processes of all users are jointly stationary and ergodic, with continuous stationary distribution.^ 
The joint fading process is slowly time-varying relative to the codeword's length, and adheres to a 
block fading channel model, which remains constant for a time block T, but is allowed to change 
in an independent identically distributed (i.i.d.) fashion from block to block. This is a valid model 
for ideally interleaved TDMA or packet-based access where each data frame "sees" an independent 
channel realization which remains constant within each frame [15, Chapter 2]. User transmissions 
to the access point are naturally frame-based, where the frame length is chosen equal to the block 
length. Having perfect knowledge of the (possibly quantized) {/ife}^i, the access point assigns 
time fractions to users via a (uplink map) message before an uplink frame. Then users transmit 
with the rate adapted to their CSI at the assigned time fractions. Let h := [/ii, . . . , hxY denote 
the joint fading state over a block. Through feedback from the access point, the K transmitters 
are assumed to know h and can vary their codewords, transmission rates and transmission times 
per block. 

A/^otation: We use boldface lower-case letters to denote column vectors and inequalities for vectors 
are defined element-wise. We let -F(h) denote the cumulative distribution function (cdf) of joint 
fading states, £'h[-] the expectation operator over fading states, f^^\x) the fcth derivative of f{x), 
(p the empty set, ^ the transposition operator, I^.j. the indicator function (I^^y = 1 if a; is true and 
zero otherwise), and := max(x, 0). 

'As with [14], our analysis can be easily extended to discrete distributions. 

3 



K 




(1) 



k=l 



III. Weighted Sum Average Rate Constraint 

We first consider the problem of minimizing total power given a weighted sum average rate 
constraint. Such a constraint may arise in a wireless sensor network, where the fusion center 
requires an aggregate rate R to perform a certain task (e.g., distributed estimation) using data 
from different users with different reward weights. Given a rate allocation policy r() and a time 
allocation policy t( ), let rfc(h) and rfc(h) denote the time fraction allocated to user k and the 
corresponding transmission rate during Tk{h). Taking into account that user k does not transmit 
over the remaining 1 — Tfe(h) fraction of time, the kth user's overall transmission rate per block is 
Tfe(h)rfe(h). If w := [wi, . . . ,wk]^ collects the rate reward weights assigned to the K users, we 
let JF^ denote the set of all possible rate and time allocation policies satisfying the average rate 
constraint Eh[J2k=i''^kTkO^)i^kO^)] > R with '^k=i^kO^) = 1> Vh. Clearly, using transmit power 
Pfe(h) during Tfe(h) fraction of time in any given block, user k can theoretically transmit with 
arbitrarily small error a number of bits/sec up to the Shannon capacity rfc(h) = B log ^1 + hh^j^^ ^ 
where B is the system bandwidth. Without loss of generality (w.l.o.g.), we assume henceforth that 
B — 1 and cr^ = 1. Again, notice that with allocated time fraction Tfe(h), the kth user's overall 
transmit power per block is Pk{ii-) — Tfe(h)pfe(h) since no power is used for 1 — Tfe(h) fraction of 
time. 

With Pk := Eh[Pk{h)], p := [Pi, • • • , Pk]^ and in accordance with the definition of the ergodic 
capacity region, we define a power region as follows. 

Definition 1: The power region for the TDMA fading channel when transmitters and the receiver 
have perfect CSI, is given by 

V{R)^ U Vtd{v{-)M-)), (2) 

(r(.),x(.))e^w 



where 



VTD{r{-),T{-)) = \p:Pk>Et, 



^fe(^) ('2^'^('^) - 1) 

hk 



l<k<K}. (3) 



If the block length is sufficiently large and the users are allowed to use different codewords 
for different fading states, it is easy to show that every p e P(-R) is feasible. Moreover, by the 
time-sharing argument, we can show that the X-dimensional power region V{R) is convex in p 
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(the proof mimics the steps in capacity region derivations [9], [11], and is omitted for brevity). 

Supposing that we assign to users different weights fx :— [//i, . . . , ^ 0, the energy-efficient 
resource allocation problem can be formulated as 

min/x^p, subject to pEV{R). (4) 
p 

Its solution p yields the optimal rate and time allocation policies, and lies on the boundary surface 
of V{R) due to its convexity. By solving (4) for all /x > 0, we can determine all the boundary 
points, and thus the entire power region V{R). When one or more of the entries of fx are zero, 
the solution to (4) corresponds to an extreme point of the boundary surface of V{R). By letting 
some of the weights approach 0, we can get arbitrarily close to these extreme points. One can 
refer to [9], [10], [11], [12], [14] for the explicit characterization of the extreme points. 



A. Full CSI and Infinite-Codebooks 

When a user can vary its codebook according to each fading state, the boundary of V{R) is 
feasible. Therefore, the problem (4) can be rewritten as 



subject to Eh Ef=i«^/kT-/fc(h)rfc(h) = R, Vh Efc=i^fc(h) = l- 



(5) 



Using the Lagrange multiplier approach, we can decompose (6) into two sub-problems. 

1) Given what we term the total rate-reward -R(h) := Y^k=i'^kTk(}^)Tk{)^) assigned to the K 
users, we determine how to distribute -R(h) among users so that the total power cost in a 
fixed state h is minimized. That is, we solve 

f J(i?(h)) := minr(h),.(h) Ef=i -".•^(2^'=('^) - 1) 
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(6) 



subject to Ef=i WfcTfe(h)rfc(h) = i?(h), Ef=i '^fc(h) = 1- 



Upon defining the rate-reward for user k as -Rfe(h) :— u'ferfc(h) and the corresponding power- 
cost~rate-reward (C~R) function 

/,(x):=^(2-/""=-l), (7) 
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we can rewrite (6) as 

J(i?(h)) := minr(h),T(h) 
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(8) 

subject to X)f=i Tifc(h)-Rfe(h) = R{h), Ylk=i M'^) = 1- 



2) Having obtained J( ) in (8), we optimize the allocation of -R(h) across the realizations of 
h, so that the total power cost averaged over all fading states is minimized; that is 



min^(h) E^[J{R{h))] - A^h[i?(h)] 

(9) 

subject to Eh[R{h)] = R 
where A denotes the associated Lagrange multiplier. 
To gain insight, we first solve (8) and (9) for two-users before generalizing to i^T-users. 
1) Two-User Case: With fi{x) and f2{x) denoting the C~R functions corresponding to users 
1 and 2, we first establish following lemma. 
Lemma 1: Supposing w.l.o.g. that Wi < W2, it holds that: 

1) ^/ ^ > 'hen Mx) < Mx), Vx > 0. 

2) If-^-< -^r, then 

f2ix) > fiix), when < x < vq, 

(10) 

f2{x) < fi{x), when x > Vq] 

where vq is the unique solution to the equation f2{x) = fi{x). 
Proof: See Appendix Al. □ 
When wi < W2, Lemma 1 asserts that if > the C~R curve fi{x) of user 1 stays 

always above f2{x). If < the two C~R curves cross each other once at vq, as shown in 

Fig. 1; hence, Vx G (0, vq), we have f2{x) > fi{x); and for x > vq, f2{x) < fi{x). Using Lemma 

1, we can characterize J(i?(h)) in (8) as follows. 

Lemma 2: For K — 2 and Wi < W2, the solution to (8) is: 

V — '^1^2' "^(-^C^)) ~ /2(-R(h)), which is achieved by the allocation Tj*(h) = 0, 
r*(h) = 0, r*(h) = 1, and r*(h) = R{\i)/w2. 
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2) If^< -^r, then 

( 

h{R{h)), if < R{h) < i?„(h); 

J{Rm = I f2{R{h)), if R{h) > R,{h); (11) 

fi{Ra{h)) + 5o(h)[i?(h) - Ra{h)], if Ra{h) <R< R,{h); 



where 



/ ,v / I / --z-z , 1 . XWihA W2-Wi\ yU2 /il 



onti so(h) '■s' solution of the equation 

f XW2h2\ 

i.e., 5f(so(h),h) = 0. The minimum cost J{R{h)) is then achieved with these policies: 

a) j/O < R{h) < i?„(h), then 

'ri*(h) = l, Tl{h) = R{h)/w^, 

T2*(h) = 0, r2*(h) = 0; 

b) //■i?(h) > Rkih), then 

'Ti*(h) = 0, rr(h) = 0, 

_r|(h) = l, r^(h) = i?(h)/«;2; 

c) ifRa{h) < R{h) < Rb{h), then 



(14) 



(15) 



(16) 



T*(h) = 1 - Ti(h), r*(h) = R,{h)/W2. 

Proof: See Appendix A2. □ 
Lemma 2 specifies the optimal power cost curve J(i?(h)) for each fading state h when ivi < W2- 
Specifically, if > J(-R(h)) is simply /2(-R(h)); otherwise, J{R{h)) comprises part of 

the /i(i?(h)) curve, the tangent line, and part of the /2(-R(h)) curve, as indicated in Fig. 1. Having 
obtained J(i?(h)), we now solve (9) to obtain the optimal resource allocation policies. 
Theorem 1: For K = 2 and Wi < W2, the optimal rate and time allocation policies with respect 
to the minimization problem (4) are as follows: 



1) If^> then 



r|(h) = l, r,*(h)= [logA*-log(^ 



(17) 



((«)lfel)/(ln2^l))"'l \ ^2-^1 
((i«2/i2)/(ln2M2))'"2 



a) if \* <i or if X* > ^ and g{X*, h) < 0, then 



log A*- log 



wihi 



Tr(h) = l, rt(h) 

T2*(h) = 0, r^(h) = 0; 

b) if X* > ^ and g{X*, h) > 0, 

Tr(h) = 0, rr(h) = 0, 

r2*(h) = l, r^(h)= [logA*-log(^ 

c) A* > ^ and g'(A*, h) = 0, then for an arbitrary Tq G [0, 1], 

Tr(h) = To*, rr(h) = [log A* - log (^) 



(18) 



(19) 



+ 



log A* — log 



In 2/12 
W2h2 



(20) 



r*(h) = l-r*, r*(h) 

Function g{x, h) 15 gjven by (13), and A*, Tq are obtained numerically by satisfying the weighted 

sum-rate constraint -£'hEfc=i '"^fe''"fe(h)r^(h)] = R. 

Proof: See Appendix B. □ 
Note that depending on A*, water-filling in (17)-(20) may result in zero transmission rate for 

the user which has been assigned the entire or part of the block. Therefore, for some fading states 

where the channel is really bad, both users should defer. Comparing (8) and (9) with [11, eqs. 

(11), (13)], we find that our power minimization yields similar "opportunistic" policies as the rate 

maximization in [11]. 

Since Tq can take any arbitrary value between and 1, the solution in (20) is not unique. 
However, under the assumption of continuous joint fading distribution density, the probability of 

g{X*, h) = is zero, and therefore case c) is an event of measure zero. Thus, after setting Tq to an 
arbitrary value in [0,1], A* can be uniquely determined by an one-dimensional, e.g., bi-sectional. 



8 



search. From Theorem 1, it is clear that in order to achieve energy efficiency over TDM A fading 
channels, most of the time we should allow one user to transmit per block. This also holds true 
in rate maximization for TDM A fading channels [11], even though the resultant time and rate 
allocation fractions are different. 

To extend our two-user results to K > 2 users, we will need definition of the convex envelope. 
Definition 2: The convex envelope f^{x) of a function f{x) is the solution to the optimization 
problem 

^'{x) = max ax + b, subject to ax + b < f{x), \/x. (21) 

a,b 

Namely, f^ix) is the boundary surface of the convex hull of the function's epigraph [16]. 

Using the definition f{x) := m.m.i<k<K fk{x), we can verify the following property: 
Proposition 1: The optimal C'^R function J(i?(h)) in (8) is the convex envelope of f{R(h.)) in 
the two-user case; i.e., J(i?(h)) = f'^{R{h)) as determined by Lemma 2 and Definition 2. 

Proof If wi < W2 and > then J(i?(h)) = /2(i?(h)) = /(i?(h)). It is trivial to 

show J(i?(h)) = /'^(i?(h)). If wi < W2 and ^ < then it is easy to show that J{R{h)) 
is convex since it has non-decreasing first derivatives for all -R(h) > 0. If the convex envelope 
/'^(/?(h)) were not given by J(i?(h)), then f'^{R{h.)) would be strictly greater than J{R{h)) for 
some R{h). Since the first and the third branches of J{R{h)) in (11) are exactly f{R(h)), for 
< R{h) < Ra{h) and R{h) > Rb{h), we must have r{R{h)) = J(i?(h)). Therefore, r{R{h)) 
can only be greater than J{R{h)) for R{h) e {Ra{h.), Rf,{h)). But since /"^(it!(h)) is convex, 
Jensen's inequality implies that its value for any Ra{h.) < R{h) < Rb(ii) can not be greater 
than the value given by the line segment connecting f{Ra{h)) with /(i?(,(h)). This leads to a 
contradiction, and thus J{R{h)) = f%R{h)). □ 

2) K-Use Case: Generalizing Proposition 1 to X > 2 users, we can show that: 
Theorem 2: For a K-user TDMA block fading channel, the optimal C'^R function J(i?(h)) at 
each fading state h is the convex envelope of f{R{h)) :— mini<fc<x /fc(-R(h)), and the optimality 
is achieved by allowing at most two users to transmit per time block. 

Proof: See Appendix C. □ 
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Although Theorem 2 asserts achievability of the optimal policies, it does not provide algorithms 
realizing these optimal allocation strategies. The latter are challenging since obtaining convex 
envelopes is generally difficult. As pointed out in [17], the rate-maximizing resource allocation for 
TD broadcast fading channels in [11] is obtained by water- filling across realizations of concave 
envelopes. Our energy-efficient resource allocation policies under a weighted sum average rate 
constraint for TDMA fading channels are obtained via water-filling across realizations of the convex 
envelopes J(i?(h)). But in order to determine J(i?(h)), we must generalize Lemma 2 to the K-usev 
case. This is accomplished with the following algorithm. 

Algorithm 1: initialization: Let u := [ui, . . . , ukq]'^ (initially set equal to the empty set 0) denote 
the set of active users during a given block, s := [si, . . . ,sko\'^ the set of slopes of (possibly 
multiple) tangent lines common to the active C~R curves (also initially set to 0), and let the 
iteration index be m = 1. 

51) Remove the C~R function of user k from J(R(h)) if 3i k such that Wk < Wi and 

> because in this case Vx > 0, fkix) < fi{x) and fk{x) will not appear in the 
expression of J(i?(h)) for the reasons we detailed in Lemma 2. Let Kr be the number of users 
remaining after such a successive elimination of their C~R functions from J{R{h)), and 
define the permutation 7r(-) such that W7r(i) < Wt^(^2) < ■ ■ ■ < WT^{Kr)- Then through successive 
pairwise comparisons and user C~R function removals, we can ensure that — — < 

52) Let the mth element of uhe Um — 7r(l). If Kr > 2, go to Step 3. If Kr < 2, all users with 
C~R functions not appearing in J{R{h)) have been removed. Then set the number of active 
users be Kq = m, Sm = oo, and stop. 

53) Forl<i<Kr-l, define 

h) := X (^».,.«, log ( ) - log [■^^^^j ) 

\ In 2 / K{i) ' 

and find the 6 satisfying g,{^,,h) = 0, log [ """"^Xt^^T ) > ° ^""^ ( TC" ) > 0- 
Let also = minj{^j}, and i* = arg minj{^j}. Remove C~R functions of users 7r{k) for 

10 



which 1 < k <i*. Increase m by 1 and return to SI). 
Using Algorithm 1, we obtain u and s. If Kq — 1, J{R{h)) is simply fui{-)', if Kq > 1, then 
J(/?(h)) comprises pieces of the curves well as the common tangent line segments 

between and I < m < Kq - 1, the slopes of which are {s„t}^'l^^ By denoting 

i?„o(h) = -R6o(h) = 0, -Ra^fjh) = -R6Ko(h) = oo, and letting -Ra^(h) and -R6^(h) be the points 
with equal first derivatives /ij„^(i?a^(h)) = fu]l+i(RbraO^)) = Sm, we can write 

fujR{h)), if i?6_,(h) < R{h) < Rajh); 

(23) 

/ujRaJh)) + Srr,[R{h) - R^Jh)], if R^Jh) <R< R,Jh). 

Arguing as in the proof of Proposition 1, it follows readily that J(i?(h)) is the convex envelope 
of f{R{h)). Once having J{R{h)), we will implement a water-filling strategy to obtain the energy- 
efficient resource allocation across realizations of h. First, for any user /c ^ u, we let T^(h) = 
and r^(h) =0. The rate and time allocation for the remaining Kq users is given as follows. 
Algorithm 2: 1) If Ko = 1, then 

'ln2/i„^' 



Am) 



log A* - log 



(24) 



2) If Kq > 1, set So = 0. Since = Sq < Si < • • • < = oo, for the given A*, there exists 
J e {1, 2, ... , Kq} such that < A* < Sj or A* = sj. 

a) If 3j such that < A* < Sj, we know J(i?(h)) = /„^.(i?(h)). In this case, we set 

'ln2 fly 



log A* — log 



(25) 



<,(h) = l, r;.(h) 

and Vi 7^ j, r:^(h) = and ^^(h) = 0. 
b) If 3j such that A* = Sj, then J(i?(h)) = (h)) + Sj[R{h) - i?a.(h)]. As in the 

two-user case, we set 

-;(h)=ro*, r;.(h)=[logA*-log(^ 



t: (h) = 1 - T* rl (h) = [log A* - log ( '"'^;^-+^ \ 
and Vi ^ J, J + 1, T;(h) = and <^(h) = 0. 
In eqs. (24)-(26), A* and Tq are chosen to satisfy the weighted sum-rate constraint 



(26) 



K 



fe=i 



R. 



(27) 
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Under the assumption of continuous joint fading distribution density, the value of Tq does not 
affect the weighted sum-rate constraint, and A* can be uniquely determined by an one-dimensional 
search, as in the two-user case. To achieve energy efficiency, we should only allow at most two 
users (and most of the time only one user) to transmit per time block in the i^-user case. Again in 
(24)-(26), water-filling may result in zero transmission rate for the user which has been assigned 
the entire or part of the time block. For some fading states, when all channels are in deep fading, 
all users should defer. Also note that Algorithms 1 and 2 are dual to those in [11] for power 
minimization under an average sum-rate constraint. 

B. Quantized CSI and Finite AMC Modes 

In this section, we provide a novel formulation and solve the energy minimization problem 
under a weighted sum average rate constraint for the finite-AMC-mode case. In practice, a user 
may not be able to support an infinite number of codebooks. Moreover, the codewords in use may 
not be capacity-achieving. It is thus worth investigating energy-efficient resource allocation for 
practical systems where each user can only support a finite number of AMC modes. Notice that 
since transmitters can transmit with a finite number of AMC modes, only quantized CSI could be 
fed back from the access point to the transmitters suffices. 

For user k G [1, K], an AMC mode corresponds to a rate-power pair {pk,hPk,i), l = 1, • • • , M^, 
where denotes the number of AMC modes. A pair {pk,uPk,i) indicates that for transmission rate 
Pk^i provided by the Zth AMC mode, the minimum received power required is pk^i- Notice that the 
minimum power pk^i may not be given by 2^'= ' — 1 as in the capacity-achieving case, and some extra 
power may be required in practice. Also, the rate pk,i is maintained with a prescribed symbol error 
probability (SEP), and Pk,i is the corresponding minimum received power under the SEP constraint. 
For this reason, we need to implicitly include the SEP constraints in our optimization. Although the 
fcth user only supports Mk AMC modes, this user can still support through time-sharing continuous 
rates up to a maximum value determined by the highest-rate AMC mode Pk,Mk- 

By setting pk^ = and pkfi = and letting -fk,i := (Pk,i - Pk,i-i) / {pk,i - Pk,i-i), we define the 
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piece-wise linear function relating transmit power with rate as 

{Pk,i-i/hk + 'yk,i{'f^k{ii) - Pk,i-i)/hk, Pk,i-i <rk{h.) < pk,i, le[l,Mk]; 
(28) 
oo, rfc(h) > pk,Mk- 

Notice that in order to support rate pk,i with channel coefficient hk, the required transmit power is 

scaled as Pk,i/hk. For practical modulation-coding schemes with e.g., M-QAM constellations and 

error-control codes, Tfe(rfe(h)) is guaranteed to be convex [2]. Using (28) to replace the power-rate 

relationship implied by Shannon's capacity formula, we can define a power region as [c.f. (2), (3)] 

V'{R)= U VUr{-),r{-)), (29) 

(r(.),T(.))e^w 

where 

^TD(r(-), = {p : Pit > [Tfe(h)T,(rfe(h))] , 1 < k < K} . (30) 

It is easy to show that the i^-dimensional V'{R) is feasible and convex. The optimization problem 
thus becomes 

min/Li^p, subject to pEV'{R). (31) 
p 

We can rewrite (31) as 

minr(.),x(.) Eh Ef=i /Xfcrfe(h)Tfc(rfc(h))] 

(32) 

subject to El, ^f^^WfeTfe(h)rjfc(h)] = ^, Vh Ef=i '^ik(h) = 1- 
As in the infinite-codebook case, we can still decompose (32) into two sub-problems. That is, we 
first optimize per fading realization, and then apply water-filling across J{R{h)) realizations to 
obtain the optimal resource allocation policies. 

Recalling the weighted rate-reward -Rfc(h) = Wkrk{h), the C~R function corresponding to the 
kth user is now /fc(i?fe(h)) := pk'^ k{Rk()^) / Wk) which is a piece-wise linear curve through the 
points {{wkpk,i-, IJ'kPk,i/hk)}fi\- Since the C~R functions are piece-wise linear, the convex envelope 
J{R{h)) of their min^ fk{R{h.)) can be obtained by simply comparing the slopes of a finite number 
of straight line segments. To this end, we implement the following algorithm: 
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Algorithm 3: initialization: Define the set of points m.k :— [{wkPk,i: IJ'kPk,i/hk), / = M^]^, 
V/c e [1, -f^']- Start with the set of slopes s := [sq, ■ ■ ■ , s^-i]^, set of rate rewards r := [i?o, ■ ■ ■ , Rm-iV, 
set of power costs c := [Co, . . . , Cm-iY with sq = Rq = Cq = 0, and let the iteration index be 
m = 1. 

51) Consider the point (xo,yo) — (-Rm-i, C'm-i)- For k — 1,...,K, if m.^ ^ 0, let the ith 
element in nife be {wkpk,rnk[i), l^kPk,mu{i)/hk)- Also, set {xk.Vk) = {wkPk,rnk{i),Pk,rnk[i)lhk)- 

52) Let Sra = miiifc: rn^,^^ and k* = arg miiifc; ^^^^ Set i?^ = Xk* and = yk*. 

a) For each k ^ k* and nifc ^ (j), k e [l)-^]> ^* = maxj {i : WkPk,rak(i) < 
tyfe*Pfc*,irn^*(i)}- If 7^ 0, then remove all {wkpk,mk(i), P-kPk,rnk{i)lhk), i e from 
rhfc, and reduce by i* . Note that in the next iteration, the original (tffcPfc.mfe(j*+i), 
PkPk,rn^{i*+i)lhk) becomes the first element {wkPk,ink{i)^ l^kqk,m^{i)/hk) of nifc. 

b) Remove {wk*Pk*,ihk*{i): l^k-Pk*,ihk*{i)/hk*) from m^* and reduce Mk* by 1. 

If iTifc — (/) (i.e., Mfc = 0), yk e [1, i^], set Kq — m and stop. Otherwise, increase m by 1 
and go to 57 ). Notice that Kq is the number of comer points of the wanted convex envelope 



An illustration example for Algorithm 3 and the resultant J(i?(h)) is shown in Fig. 2, where 
each of the two users can support three AMC modes. In SI of Algorithm 3, we first compare the 
slopes and //2P2,i/ (^21^2^2,1) and find /i2P2,i/(^2«^2P2,i) < A«iPi,i/(^iiyiPi,i)- 

Therefore, we let Ri — ^2^2,1 and Ci = //2P2,i/^2- Since both wipi^i and wipi^2 are less than 
W2p2,i, we remove the first two AMC modes of user 1 from rrii in S2 a); whereas we remove the first 
AMC mode of user 2 from m2 in S2 b). In SI of the next iteration, we compare the slope between 
points (wipi.s, P'iPi,z/hi) and (^2^2,1, A«2P2,i//i2), with the slope between points (^2^2,2, A«2P2,2//i2) 
and (1^2^2,1, A*2P2,i/^2)- In this case, we should set R2 — wipi^s and C2 = A*iPi,3/^i- In S2, we 
remove (w2P2,2, P2P2,2/h2) from m2, and remove (wipi.3, jiipi^^jhi) from rni. In the last iteration, 
we obtain R^ — ^2^2,3 and C3 = P2P2,3/h2, and J(i?(h)) is determined. 





(33) 
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Having determined J{R{h)) as in (33), we implement water-filling across realizations of J(i?(h)) 
to derive the optimal resource allocation policies. Different from the infinite-codebook case where 
all J(i?(h)) have continuous slopes, here J{R{h)) are piecewise-linear. Therefore, water-filling 
should take into account the finite number of slopes of J{R{h)). A somewhat related problem 
was dealt with in [12], where water-filling over some piecewise-linear concave functions was used 
to determine the boundary surface of the outage probability region. But energy-efficient resource 
allocation policies for finite- AMC-modes were not considered in [11] and [12]. 

Recall that Kq and all entries of s, r and c are functions of h. Since every point of the convex 
envelope J{R{h)) can be achieved by time-sharing between points (it!^(h), C^(h)), finding the 
optimal resource allocation strategies is equivalent to solving the following minimization problem: 



min^(h) El, Em^i '^'»(h)Cm(h) 



subject to El, E^Li fmih)RUh) > R, ^2^1 ^m(h) = 1. 



(34) 



^^s^o 



Theorem 3: If the optimization (34) is feasible, by setting sxo+iih) = oo, then Vh, we have the 
optimal solution {T^^(h)}^!^^, and thus the optimal allocation policies T^(h) and r^(h) (1 < k < 
K) for the original problem (31) as follows. • 

1) If \* < si(h), then f^(h) = 0, Vm = 1, • • • , Kq; consequently, T^(h) = and r^(h) = 0, 
Vk = l,...,K. 

2) If 3m* e {1, 2, . . . , Kq} so that Sm*(h) < A* < Sm*+i{h.), then f^*{h) = 1, and f^{h) = 0, 
Vm m*, 1 <m < Kq. This implies that if (i?^*(h), Cm*{}oL)) belongs to user k*, then 

T,^(h) = l, rl.{h)^Rm^{h)/wk*; (35) 

and T^(\i) = and r*(h) = 0, ^ k*, 1 < k < K. 

3) If 3m* e {1, 2, ... , Kq} so that X* = s„*(h), then f^*(h) = Tq, f^*_i(h) = 1 - Tq, and 
f^(h) = 0, Vm 7^ m*,m* — 1, 1 < m < Kq. This implies that if (i?^*(h), C^*(h)) and 
{Rm*-i{h),Cm''-i{ii)) belongs to users i* and j*, respectively, then 

r^(h)=To*, r*.(h) = i?^*(h)M*, 

(36) 

7;.(h) = 1 - To*, (h) = R^._,(^h)/w,,; 
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and T^{h.) = and r^(h) = 0, V/c 7^ 1 < k < K. Note that if m* = 1, we let user i* 

transmit with Tq fraction of time and leave the channel idle for the remaining 1 — Tq fraction 
of time. If i* = j*, then we let the same user transmit with two AMC modes, one for Tq 
fraction and the other for 1 — Tq fraction of time. 
In eqs. (35) and (36), Tq should satisfy the average rate constraint 



The results of Theorem 3 are analogous in form with those in [12, Theorem 3], which is a 
generalization of [13, Lemma 3]. But note that the latter deal with the delay-limited/outage capacity; 
while our results are for energy-efficient resource allocation using finite-AMC-modes, a subject not 
considered in [12], [13]. In Theorem 3, A* is the water-filling level. For energy efficiency, we should 
let the first derivatives J^^\R{h)) = A*. However, since J(i?(h)) entails a finite number of slopes, 
equality can not be always achieved. Thus, our strategy is to select the largest J^^\R{h.)) < A*. 
When the largest J^^\R{h.)) — s^*(h) < A*, i.e., s^*(h) < A* < s^«+i(h), since the user(s) 
transmit more efficiently than the required power level, we should allow transmission(s) with peak 
rate given by i?m*(h)/wfc*. When A* = Sm.(h), users transmit as efficiently as required; thus 
arbitrary time division suffices. When A* < Si(h), no transmission can be carried out as efficiently 
as required, and all users defer during this fading state. As in the infinite-codebook case, we 
should allow at most two users to transmit per time block. However, here we also allow one user 
to transmit in a time-sharing fashion with two AMC modes during some time blocks. 

C. Discrete-Time Allocation 

So far we have derived energy-efficient resource allocation strategies under the assumption that 
the time fraction assigned to each user can be any real number in [0,1]. In practical TDM A 
systems, time is usually divided with granularity of one time unit (slot), which is determined by 
the available bandwidth. Therefore, the transmission time assigned to each user per time block 
has to be an integer multiple of a "slot". Take the well-known GSM system as an example [15, 
Chapter 4]. In each narrowband channel of 200 KHz, time is divided into slots of length 577 //s 




(37) 



Proof: See Appendix D. 



□ 
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and each uplink frame is shared by the users in a time-division manner, and consists of 8 slots. 
Upon regarding an uplink frame as a frequency flat-fading time block, we can thus only assign 
each user a time fraction which is an integer multiple of 1/8. However, this extra discrete-time 
allocation constraint does not affect the derived energy-efficient allocation policies. In our policies, 
most of the time we should assign the entire time block to a single user. When we occasionally 
allow two users to transmit (or allow one user to transmit with two AMC modes), the time division 
between them can be arbitrary. Therefore, our energy-efficient policies can be easily adopted by 
practical systems which only allow discrete-time allocation among users. 

IV. Average Individual-Rate Constraints 

In resource allocation for multi-access channels, other than a weighted average sum-rate con- 
straint, a more general setting is when each user has an individual average rate requirement. We 
next consider power minimization under such individual rate constraints. Without the rate-reward 
weight vector w, we let T' denote the set of all possible rate and time allocation policies satisfying 
the individual rate constraints {-Eh[TA;(h)rfc(h)] > -R^}^^ and ^Y^=\ ^k^) = 1, Vh. Upon defining 
f := [^1, . . . , Rk]^, the power region under the individual rate constraints is [c.f. (2)] 

p(r)= U VrD{r{-),T{-)), (38) 
{r(.),x{.))e^' 

where Vtd{y{ )-, '''{■)) is defined as in (3). Again, if the block length is sufficiently large and the 
users are allowed to use an infinite number of codebooks, it is easy to show that every point 
in V{v) is feasible and the region is convex. With power cost weights /* := [//i, . . . ^HkY, the 
energy-efficient resource allocation policies solve the optimization problem 

min/Li^p, subject to p e 'P(f). (39) 
p 

The solution p yielding the optimal rate and time allocation is on the boundary surface of V{y) due 
to its convexity. By solving (39) for all /a > 0, we determine all the boundary points, and thus the 
whole power region ■P(f). Again, we will explicitly characterize the optimal resource allocation 
policies and the resultant boundary point for /it > 0. By letting some of the power cost weights 
approach 0, we can get arbitrarily close to the extreme points. 
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(40) 



A. Infinite-Codebooks 

When the user can utiUze an infinite number of codebooks to achieve channel capacity per fading 
state, the optimization in (39) is equivalent to 

minr(.),.(.)£;h[Ef=i/^.^(2-('^)-l) 
subject to Eh [Tfc(h)rfc(h)] > Rk, k = 1, 

Ef=iT.(h) = l, Vh. 
The counterpart of Theorem 2 in this case is given by: 

Theorem 4: For any /J, > 0, there exists a A* := [A^, . . . , A^]''^ > 0, and optimal rate and time 
allocation policies r*(-) and t*(-) in (40), such that for each h, r*(h) and T*(h) solve 

niinr(h),r(h) Ef=i^fe(h)/fc(rfe(h)) 

(41) 

subject to Y^k=i ^kO^) = 1' 
where now /fc(rfe(h)) := |^(2'''='^^) ~ 1) ~ '^fc^fc(h). Since /fc(rfe(h)) is convex in rk{h), it attains its 
minimum at rfe^inin(h) := log A)^ — log . Moreover, we have r*(h) and T*(h) as follows. 

1) If functions fk{rk,minO^)), k — 1, . . . , K, have a single minimum /j(rj,inin(h)), i.e., 

i = arg min fk{rk,mm{^))^ 

k 

then yk^i, ke [1, K], r^(h) = 0, T;{h) = 0, and 

r*(h)=n,^in(h), T*(h) = l. (42) 

2) If functions fk{rk,minO^)), k = 1,...,K, have multiple minima {/i^ (rj^,min(h))}^^i. then 

r*.(h)=r.^,„..n(h), T-.(h) = r;, (43) 

arbitrary Ej=i "^Z = 1- h' k K], r^(h) = and T^ih) — 0. 

In (42) and (43), X* and {t*}^= I are obtained by satisfying the individual rate constraints 

E,,[T;{h)rl{h)]^Rk, k^l,...,K. (44) 
Proof: See Appendix E. □ 
If we regard /^(rfc inin(h)) as a channel quality indicator (the smaller the better) for user k. 
Theorem 4 asserts that for each time block, we should only allow the user with the "best" channel 
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to transmit. When there are multiple users with "best" channels, arbitrary time division among them 
suffices. Therefore, our resource allocation strategies are "greedy" ones. Note that fk{rk,mmO^)) 
contains A^. This implies that the user having smallest fk{rk,mm{ii)) actually has the rate-constraint- 
controlled "best" channel. 

Rate maximization under individual power constraints was pursued in [9], where it was shown 
that superposition codes and successive decoding should be employed and that greedy water- 
filling based on a polymatroid structure provides the optimal resource allocation. In our power- 
efficient TDMA setting, we do not have such a polymatroid structure. Albeit in different forms than 
those in [9], our strategies can be also implemented through a greedy water-filling approach. To 
obtain the optimal allocation policies in Theorem 4, we need to calculate the Lagrange multiplier 
vector A*. Although a i^T-dimensional search can be used to directly compute A* from (44), it is 
computationally inefficient when K is large. Next, we show that an iterative algorithm from [9] 
can be adopted to calculate A*. Before that, by the strict convexity of exponential functions and 
the fact that non-uniqueness of the time allocation t*( ) occurs with probability when -F(h) is 
continuous, we can argue as in [9, Lemma 3.15] to establish the following lemma: 
Lemma 3: Given a positive power weight vector /x, there exists a unique p* which minimizes 
fM^p, and there is a unique Lagrange vector A* such that the optimal rate and time allocation 
satisfy the average individual rate constraints. 

To gain more insight, let us look at a special case where the fading processes of users are 
independent. If ) stands for the cumulative distribution function (cdf) of user k's fading channel. 



Corollary 1: If the fading processes of users are independent, the optimal solution p* to (39) for 
a given fi > can be obtained as 



define gi{x) : 
gk{z) = 0. 



In 2 



X 



K log i^j-:, and let Si^kiz) denote the solution to gi,k{x, z) :— gi{x) 




(45) 



where the vector A* is the unique solution to the equations 





(46) 
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Proof: By definition, we have 



.(h) = 



log XI - log 



In 




0<hk< ^ 



Kh 



±k. _ M _ \* Ino- —J 
In 2 hk 6 ln2/ifc 

0, 



0<hk< ^ 



and 



^/fc(rfc,min(h)) _ J i (K 1112) ' > 



0, 



0<hk< ^ 



Since 



9 f kirk ,min (h)) 
dhk 



< 0, we have from (49) [c.f. Theorem 4] 



^h[r,*(h)r*(h)] = i?h[r*(h)I 



poo 

/ log 



{/i(''i,min(h))>/fc(rfe_min(h)), Vi} 



and 



t2'"k,mUh) _ l)I|^,(^^^^.^(h))>/fe(r-(h)), Vi} 



In 



(47) 



(48) 



(49) 



(50) 



(51) 



□ 



The corollary thus follows from (50) and (51). 

What left to obtain the optimum {Pjf}^]^, is to specify {Xlj^^i- We accomplish this by 
modifying the corresponding iterative algorithm in [9]. 

Algorithm 4: Let A(0) > be an arbitrary initial rate-reward vector. Given the Zth iterate 

the (Z + l)st iterate + is calculated as follows. For each k e [1, -f^], Afe(Z + 1) is the unique 

rate-reward weight for the fcth user such that the average rate of user k is Rk under the optimal 
rate and time allocation policies, when the rate-reward weights of other users remain fixed at A(/). 

In the independent fading case, Afe(Z-l-l) is the unique root XI of (46), which can be numerically 
solved if the fading statistics are known. Let f (A(/)) := [Ri{X{l)), . . . , Rk{X{1))]'^ and p(A(Z)) := 
[Pi(A(/)), . . . , Pk{X{1))]^ denote the rates and powers of users given by Theorem 4, respectively, 
when the Lagrange multiplier is X{1). We can then prove that: 
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Theorem 5: Given the average rate constraint r, ifp* is the optimal power vector corresponding 
to the cost vector /x, and A* is the rate-reward vector satisfying (44), then 

A*, as / I oo; (52) 

and hence p(A(/)) — ^ p* and r(A(/)) — > f. 

Proof: See Appendix F. □ 
Algorithm 4 is also widely used for power control in IS -95 CDMA systems and its convergence 
is guaranteed [15, Chapter 4]. For the power minimization in our TDM A setting. Theorem 5 
ensures convergence of Algorithm 4 in finding the optimal Lagrange multiplier A*. The proof is 
analogous to that of [9, Theorem 4.3], except for some necessary modifications. With Theorem 
4 and Algorithm 4, we determine the energy-efficient rate and time allocation policies under 
individual rate constraints. Our policies are greedy in the sense that most of the time we allow a 
single user with the "best" channel to transmit, and occasionally we assign time to multiple users 
with "equally best" channels at a given time block. Note that water-filling may result in no user 
transmissions for some fading states, where all channels are in deep fading. 

B. Finite AMC Modes 

Next, we investigate optimal resource allocation under individual average rate constraints for 
the case when each user can only support a finite number of AMC modes. As in Sec. ni-B, for 
user k G [1, K], an AMC mode corresponds to a rate-power pair {{pk,i,Pk,i)}iL\, where Mk is the 
number of AMC modes. By time-sharing, a user can support continuous rates up to a maximum 
value determined by the highest-rate AMC mode Pk,Mk- By defining Tk{x) as in (28), the new 
power region is given by 

U ^TD(r(-),r(-)), (53) 
(r(.),x(.))e^' 

where J^' denotes the set of all possible rate and time allocation policies satisfying the individual 

rate constraints, and P-^dI^I')? '''(■)) is defined as in (30). It is easy to show that the region V'{f) 
is feasible and convex. The optimization problem thus becomes 

min/ii^p, subject to p G ^'(r). (54) 
p 
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Using Tk{x), problem (54) is equivalent to 



minr(.),-r(.) Eh 



Ef=i/^fc'rfc(h)Tfc(rfe(h)) 



subject to Eh [Tfe(h)rjfc(h)] > Rk, k^l,...,K, (55) 

Ef=ir.(h) = i, vh. 

Since every point of Tfc(rfe(h)) can be achieved by time-sharing between points {pk,hPk,i/hk), 
finding the optimal resource allocation strategies for (55) is equivalent to solving 



nimT.(h) 2^fe=i Eh 2^1=0 l^k 



subject to Eh Ylfio^k,i(ti)Pk 



> i?fc, A; = 1, 



(56) 



Ef=i ES^M(h) = l, Vh. 
The counterpart of Theorem 3 under individual rate constraints is now: 

Theorem 6:Iff is feasible Vh, we have the optimal solution T^^(h) (k e [Ij-ft'], I € [0, M^]) to 
(56), and subsequently the optimal allocation r^(h) and T^ih) for (55) as follows. Given a positive 
X* :— [A*, . . . , A^]-^, for each fading state h, let II :— max {I : iJ>klk,i/hk < A^} (II — if no 
such I) and Ck,i := likPk,i/hk, and define </?fc(h) := Ck,ii - KPk,ii- 

1) If {ipk{}i)}k=i have a single minimum <^i(h), i.e., i = arg miiifc (^fc(h), then fip = 1 and 



all other fk,i — 0. Consequently, 



r*(h)=p,., r;(h) = l; 



(57) 



and ^ky^i, ke [1, K], r^(h) = and T^{h) = 0. 

2) /jf {v^/c(h)}^]^ /zave multiple minima \^Lpi.(\i)y^.^^, then fi.^i* = t* with arbitrary Ej=i '^j ~ 
1, and all other fk,i — 0. Consequently, 



<(h)=r;, 



(58) 



anJ VA; ij, k e [1, K], r^(h) = and T^{h) = 0. 
In (57) and (58), X* and {t*}^^^ should satisfy the individual rate constraints 

Eh[T;{h)rl{h)]^Rk, k^l,...,K. 

Moreover, X* is almost surely unique and can be iteratively computed by Algorithm 4. 



(59) 
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Proof: See Appendix G. □ 
Theorem 6 shows that our policies with finite number of AMC-modes are still greedy ones. 
For user k at fading state h, the C~R function is given by fk{rk(h.)) := /ifcTfc(rfc(h)) — A^,rfc(h). 
It is clear that this function attains its minimum at pk,i*. Then comparing the channel quality 
indicators (pkO^) for all the users, we determine which users have the "best" channel and assign 
resources accordingly. Note that when II — 0, the user should remain silent. If at a fading state, 
this user happens to have the "best" channel, we should let all users defer at this block. When 
i^klkrjhk = K' the whole line between (pfe,;*-i,Pik,/*-i//iik) and {pk,ii^Pk,iilhk) in Tfc(rfc(h)) 
achieves the minimum of /fe(rjk(h)). Although in Theorem 6 we let user k transmit (if permitted 
by the policies) with rate pk^ip the complete solutions should allow this user to transmit with 
arbitrary time-sharing between Pk,ii-i and Pk,ip as the optimization under a weighted sum-rate 
constraint in Theorem 3. Summarizing, our greedy polices may result in no transmissions, user(s) 
transmitting with one AMC mode, or user(s) transmitting in a time-sharing fashion with two AMC 
modes per fading block. 

For the special case where the fading processes of users are independent, let Si^k{z) denote the 
solution to g'i^ki^, z) :— (pi{x) — ipk{z) — 0. Note that (pi{x) is also a function of A*. Using Theorem 
6 and mimicking the proof of Corollary 1, we can establish the following corollary. 
Corollary 2: If the fading processes of users are independent, the optimal solution p* to (54) for 
a given fi > can be obtained as 

POO 

= ^^^''^/^^ n F^{s^,k{^)) dFk{z), k^l,...,K; (60) 

where the vector A* is the unique solution to the equations 



f 



Pk,ii\\Fi{si,k{z))dFk{z)^Rk, k^l,...,K. (61) 
Some comments are now in order: 1) in the finite- AMC-mode case, the rate pk,i of user k is 



maintained with a prescribed SEP, and Pk,i is the corresponding minimum received power under 
the SEP constraint; 2) our policies for the finite-AMC-mode case could only require quantized 
CSI at the transmitters; and 3) following the arguments of Sec. III-C, it is clear that the derived 
policies for infinite codebooks and finite number of AMC-modes apply to systems where only 
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discrete-time allocations are allowed among users. 

V. Frequency Selective Channels 

The optimal resource allocation policies in the previous sections are derived for frequency- 
flat block fading channels encountered with narrow-band communications. In this section we 
extend our results to frequency-selective fading channels, which are often encountered in wide-band 
communication systems. 

Supposing that the channel varies very slowly relative to the multipath delay spread, it can be 
decomposed into a set of parallel time-invariant Gaussian multi-access channels in the spectral 
domain [18]. We consider an K-user spectral Gaussian block fading TDM A channel with continu- 
ous fading spectra uj), i?2(/, w), . . ., Hxif, uj), where frequency / ranges over the system 
bandwidth and a; is the fading state at a given time block. Let rfc(/, w) and ta;(/, denote the 
rate and fraction of time allocated to user k at frequency / and fading state u). For the weighted 
sum-rate constraint optimization, the constraint is now given by 

K 

dF{uj)>R, J]Tfc(/,u;) = l, Wf,uj. (62) 

fe=i 

Let the set .Fw consist of all possible rate and time allocation policies satisfying (62), and take the 
infinite-codebook case for illustration. The power region for this TDMA channel is given by 

U \p:Pk> [\f {2^'^^''^^ - 1) df\ dFicc), l<k<K]. (63) 

For a finite number of AMC-modes and the individual rate constraint optimization, we can similarly 
define the corresponding power regions. Subsequently, the optimal resource allocation strategies 
can be obtained from the previous results by replacing the fading state h with the frequency and 
fading state pair (/, cj) to determine power regions for frequency-selective channels. That is, we 
should employ the previous allocation policies for each (/, w), and then implement water-filling 
across both frequency / and fading state u? realizations to determine A* (or vector A*). 

VL Numerical Results 

In this section, we present numerical results of our energy-efficient resource allocation for a 
two-user Rayleigh flat-fading TDMA channel. The available system bandwidth is B — 100 KHz, 
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K 

.k=i Jf 



and the AWGN has two-sided power spectral density A^o Watts/Hz. The user fading processes are 
independent and the state h^, k — 1,2, is subject to Rayleigh fading with mean hk. Clearly, the 
signal-to-noise ratio (SNR) for user k is defined as hk/{NoB). 

Supposing hk/{NoB) = dBW, /c = 1,2, we test energy-efficient resource allocation under 
a weighted sum average rate constraint R — 200 Kbits/sec, for two different sets of rate-reward 
weights: i) wi — 1, W2 — 1, and ii) wi — 1, W2 — 2; and the resource allocation under two 
different sets of individual rate constraints: iii) -Ri = 100 Kbits/sec, R2 = 100 Kbits/sec, and iv) 
^1 = 100 Kbits/sec, R2 = 50 Kbits/sec. Fig. 3 depicts the power regions of the Rayleigh fading 
TDMA fading channels for the infinite-codebook case. It is seen that power regions I and HI under 
the weighted sum rate constraint i) and under individual rate constraints iii) are symmetric with 
respect to the line P2 = Pi. Since the individual rate constraints can be seen as a realization 
of the weighted sum-rate constraint, i.e., wiRi + W2R2 — R, the power region I contains power 
region IE. It is clear that when fii — 112, due to the symmetry in channel quality and rate-reward 
weights between the two users, the optimal resource allocation should result in ^1 = R2 under 
the weighted sum rate constraint. For this reason, the two power regions touch each other in this 
case. The relation between power regions II and IV under the weighted sum average rate constraint 
ii) and under individual rate constraints iv), are similar. They are not symmetric with respect to 
P2 = Pi due to the unbalanced rate-reward weights or individual rate constraints. Power region II 
contains power region IV, and the two regions touch each other at one point. 

For the finite-AMC-mode case, we assume henceforth that each user supports three M-ary 
quadrature amplitude modulation (QAM) modes: 4-QAM, 16-QAM and 64-QAM. For these rect- 
angular signal constellations, the SEP is given by [19, Chapter 5] 



where = 2 (l - ^) g (^^^0^ , and Q{x) := J^{l/V2^)e-y'/'dy is the Marcum's 



Q-function. From (64), we determine the rate-power pair {{pk,h Pk,i)}f=i for user k — 1,2. The 
corresponding power regions I-IV under the constraints i)-iv) for this finite-AMC-mode case with 
prescribed SEP = 10^^ are shown in Fig. 4. Similar trends as in Fig. 3 are observed. However, the 
power regions shrink since more power is required to achieve the same transmission rate in the 
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(64) 



finite-AMC-mode case relative to that in the infinite-codebook case. 

Supposing hi/{NoB) — 10 dBW and hi/{NoB) — dBW, we also test our energy-efficient 
resource allocation under the same four sets of rate constraints i)-iv). The power regions for the 
infinite-codebook case and the finite-AMC-mode case are shown in Fig. 5 and Fig. 6, respectively. 
Since the first user has a significantly better channel (i.e., higher SNR) than user 2, the required 
transmit power of user 1 is much lower than that of user 2 most of the time. Except this, the results 
are similar to those in Figs. 3 and 4. 

We next compare the derived energy-efficient resource allocation with two alternative resource 
allocation policies. Policy A assigns equal time fractions to the two users per block. Then each user 
implements water-filling separately to adapt its transmission rate at each assigned time fraction. In 
policy B, each user is assigned equal time fraction and transmits with equal power per block. Fig. 7 
depicts the power savings of our optimal policies under two different sets of rate constraints i) and 
iii), over the policies A and B for the infinite-codebook case when two users have identical SNRs. 
It is seen that when the ratio of two users' power cost weights is far away 1, our optimal policies 
under a weighted sum-rate constraint can result in huge power savings (near 20 dB) over the other 
two sub-optimal polices. However, in this case the optimal policies under individual rate constraints 
only have a small advantage (around 3 dB) in power savings, over the sub-optimal policies. This 
is because with the weighted sum average rate constraint, we can employ more flexible policies 
in time and rate allocations. From Fig. 7, we also observe that the separate water-filling in policy 
A only achieves small power savings (less than 1 dB) over the equal power strategy in policy B. 
Fig. 8 depicts the same comparison for the finite-AMC-mode case. The same trends are observed. 
However, in this case separate water-filling in policy A achieves considerable power savings (4 
dB) over the equal power strategy in policy B. Fig. 9 depicts similar power savings for the infinite- 
codebook case under two different sets of rate constraints ii) and iv), when two users have 10 
dB in SNR difference. Similar observations are obtained. But note that the optimal policies under 
individual rate constraints can also achieve large power savings (near 9 dB), over the sub-optimal 
policies. In a nutsell, our energy-efficient resource allocation policies may indeed result in large 
power savings. 



26 



VII. Concluding Remarks 

Given full or quantized CSI at the transmitters, we derived energy-efficient resource alloca- 
tion strategies for TDM A fading channels. For energy minimization under a weighted average 
sum-rate constraint, the optimal allocation policies are given by water-filling over realizations of 
convex envelopes; whereas for energy minimization under average individual rate constraints, the 
optimal strategies perform greedy water-filling. Comparing these two strategies for two different 
optimizations, we find that the first approach requires one to characterize the convex envelope of 
the minima of C~R functions per fading state, but the associated scalar Lagrange multiplier A* can 
be easily obtained by one-dimensional search. Greedy water-filling simply computes and compares 
the channel quality indicator functions of individual users and then picks the user(s) with best 
channel(s) to transmit per block; however, we need to iteratively compute the associated vector 
Lagrange multiplier A* (by Algorithm 4). 

An interesting feature of our energy-efficient resource allocation strategies should be stressed. 
According to our policies, we can let the access point (which naturally has full CSI) decide the time 
allocation and feed it back to users via uplink map messages, as in e.g., IEEE 802.16 systems [20]. 
Then given the Lagrange constant A* (with the weighted sum-rate constraint) or vector A* (with 
the individual rate constraints), the user only needs its own CSI to determine the transmission 
rate at the assigned time fraction. If uplink and downlink transmissions are operated in a time- 
division duplexing (TDD) mode, the users can even obtain their own CSI without feedback from 
the access point. Together with the fact that the access point needs only a few bits to indicate the 
time allocation (since most of the time we should allow only one user to transmit), this feature is 
attractive from a practical implementation viewpoint. 

As far as future work, it is interesting to study energy minimization over fading channels with 
delay constraint and/or using quantized CSI throughout. Our energy minimization for finite-AMC 
modes only requires quantized CSI feedback from the access point. And delay-constrained energy 
minimization may be seen as the dual problem of the delay-limited capacity maximization in [10]. 
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Extensions to these two directions are currently under investigation.^ 

Appendices 

A. Proof of Lemmas 1 and 2 

Al. Proof of Lemma 1: The kth. derivatives of f2{x) and fi{x) are given by 



= ,,ku , /{ ^(x) = —3-— 2-1 . (65) 



1^ ^Ji^ - ^ik' ^^^^^ < 2^ for wi < W2, we have /i^^(x) < fi^\x), Vx > 0. Since 

/2(0) = /i(0), we readily infer that f2{x) < fi{x), Vx > 0. 
2) If ^ < since wi < W2, there exists k > 1 such that ^^^^^ < Let kg = 

minfc argj^i^^ < Together with 2^ < 2^, we have /f "^(x) < /f"^(a;),Va: > 

0. Since /^-^)(0) > fi''-'\o) and/^)(x) < we infer that 

starts positive and with a negative slope it crosses the x-axis at some point Vko-i, hence, 

ft"~'\^) > fi'°~'\x) when 0<x< Vk,-i, 

(66) 

< fi'°-'\x) when X > Vko-i; 
Using (66) and /i^"~^^(0) > /f ""^^0), we obtain similar results for fi''°~^\x) and /f ""^^(a;) 
with Vko-2 > Vko-i. By induction, we therefore deduce that 

f2ix) > fi(x) when < x < vq, 

(67) 

f2{x) < fi{x) when x > t;o; 
where vq is the unique solution of the equation f2{x) — fi{x). 
A2. Proof of Lemma 2: For notational brevity, we drop the dependence of R, Rk, Tk (k = 1, 2), 
Ra, Rb, and sq on h. And we let ti — r and T2 — 1 — t. 

1) When wi < W2 and > we have from Lemma 1 f2{x) < fi{x), Vx > 0. We wish 
to solve (8) under the constraint tRi + {1 — t)R2 — R with -Ri > and R2 > 0. Then the 
cost function satisfies 

T/i(i?i) + (1 - T)/2(i?2) > Tf2iRl) + (1 - T)/2(i?2) > /2(i?), (68) 

^The views and conclusions contained in this document are those of the authors and should not be interpreted as representing 
the official policies, either expressed or implied, of the Army Research Laboratory or the U. S. Government. 
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where the last inequality is due to the convexity of f2{R), and the equalities are achieved 
when T — 0, Ri — and R2 — R. Inequality (68) clearly shows that the minimum in (8) is 
achieved; i.e., J{R) = f2{R)- 
2) When wi < W2 and -^J^ < fi{R) and f2{R) intersect as depicted in Fig. 1. Similar 
to [11, Proof of Lemma 1], we can specify two points Ra and Rj, such that fi^\Ra) — 
f2'\Rb) — So; i.e., at these two points first-order derivatives of fi(R) and f2(R) are equal: 

/■•"(fi.) = i^^2* = /<'>TO = !!i^2S = ., (69) 

Wihi W2h2 

The Ra, Rb expressions are obtained after solving (69) for Ra and R^. Since Sq is the slope 
of the common tangent line of the curves fi{R) and f2{R), we also have 

^ _ /2 (/?/.) -A 
-Kb — rCa 

Solving (69) and (70), we obtain sq as the solution to g{x, h) = 0, where g{x, h) is given 
by (13), and Ra and i?5 are as in (12). If < i? < 7?^ or i? > Rf,, then J{R) simply 
equals to fi{R) or f2{R)- If Ra < R < Rb, then J{R) should take the values between 
fi{Ra) and f2{Rb) on the straight line, and J{R) can be achieved by time-sharing; i.e., 
R = rRa + (1 — T)Rb, and J{R) = Tfi{Ra) + (1 — T)f2{Rb)- The optimal resource allocation 
per fading state is thus given by (14)-(16). 

B. Proof of Theorem 1 

Let i?*(h) denote the optimal total rate reward assigned to fading state h, and r^(h) and r^(h) 
the corresponding optimal rate and time fraction allocated to user A; = 1, 2. To solve (9), we rely on 
the Karush-Kuhn-Tucker (KKT) conditions [16, Chapter 5]. Let J^^)(i?(h)) denote first derivative 
of J(i?(h)) with respect to i?(h). Taking the partial derivative of E^[J{R{\v))] - XE^[R{h)] in (9) 
with respect to -R(h) for a fixed h, only J^^\R{h)) — A survives since all -R(h') for realizations 
h' h are regarded as constants. We thus have at the optimum that if R* (h) > 0, 

J«(i?*(h)) = A*. (71) 
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fl(h) 



1) If ^ > ^> then J{R{h)) - f2{R{h)) from Lemma 2, and thus J^^\R{h)) = 
Substituting the latter into (71) and recalling that -R(h) > 0, we find 



R*{h)=w2 



log 



\*W2h2 

In2/i2 



W2 



+ 



log A* — log 



In 2yU2 

W2h2 



(72) 



which yields the optimal rate and time allocation in (17). 



2) If 



Ml 



< 



At2 



wjhi W2h' 

follows from Lemma 2 that 



and we recall that sq is the slope of the straight line segment in Fig. 1, it 



if < R*{h) < Ra, A* = jW(i?*(h)) < So; 
if Ra < R*{h) < Rb, A* = J(^)(i?*(h)) = So; 
if R*(h) > Rb, X* = J«(i?*(h)) > So; 



(73) 



where Ra, Rb and sq are functions of h, given in Lemma 2. Now using (73), we arrive at: 



+ 



a) If A* < So, then J{R{h)) = /i(i?(h)), and thus R*{h) = w, [logA* -log^ 
which in turn yields the optimal rate and time allocation in (18). 

b) If A* > So, then J(R(h)) = f2(R(h)), and thus R*(h) = W2 [log A* - log ^ 
which in turn yields the optimal rate and time allocation in (19). 

c) If A* = So, then J{R{h)) = fi{Ra) + so{R{h) - Ra), and R*{h) can be any value 
between Ra and Rb. In this case, we obtain the optimal rate and time allocation (20). 

Since Ra < Rb, we have from (12) 



. SoWihi , SoW2h2 

Wi log -j — < W2 log ■ 



Wi log 



xw-i h 



In 2/^1 In 2/^2 

Moreover, since g'(so,h) = 0, and g^^\x,h.) = W2log^|^ 

(74) that 

5f(^)(x, h) < 0, for < X < ^; 

g^^\x,h) > 0, for x > ^. 
Therefore, we have: 

a') If A* < ^ or A* > ^ and ^^(A*, h) < 0, then A* < sq. 
b') If A* > C and g{X*, h) > 0, then A* > so. 
c') If A* > ,^ and g{\*, h) = 0, then A* = so. 
Combining a'), b'), c') with a), b), c) proves the second part of Theorem 1. 



(74) 
it follows from 

(75) 
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C. Proof of Theorem 2 

Since f{x) :— mmi<k<K fk{x), it is easy to show that 

K K 

J(i?(h)):= mill Vr,(h)/,(i?,(h))> min Vr,(h)/(it:,(h)). (76) 

r(h),T(h) ^ r(h),T(h) ^ 

fe=l fe=l 

By the definition of the convex envelope f'^{x), we have 

K K f ^ \ 

Y,rk{^)f{Rk{^)) > $]rfc(h)f(i?,(h)) > f $]rfc(h)i?,(h) := f(i?(h)); (77) 

k=l k=l \k=l / 

where the last inequality is due to the convexity of f'^{x). From (76) and (77), it is clear that 

J{R{h)) > f(i?(h)). 

On the other hand, according to the definition of f'^{x), for any given -R(h) there exists -Ra(h) 
and -R;,(h) (possibly -Ra(h) = -^^(h)) and r* such that 

r*i?„(h) + (1 - T*)R,{\i) = i?(h); T*f{Ra{\v)) + (1 - T*)f{R,{\i)) = f(i?(h)). (78) 

Eq. (78) shows that we can achieve the equality J{R{h)) = /''(i?(h)) by assigning 

T*(h) = T*, r*(h) = it:„(h)M 

(79) 

T;(h) = l-T*, r*(h) = i?,(h)/w,- 
to users i and j which satisfy /i(i?a(h)) = /(i?a(h)) and fj{Rb{h)) — /(i?{,(h)). This completes 
the proof. 

D. Proof of Theorem 3 

Letting M :— Yl^=i ^k-, the value of Kq per fading state may take any integer in [1, M]. We 
define Hko as the set of all fading states h for which Algorithm 3 yields Kq e [1, M]. We further 

set sxo+i(h) = cxD, VA > 0, and define for m = 0, 1, . . . , Kq, the sets 

HmiK^, A) := {h : h e Hk,, s^(h) < A < s^+i(h)} , (80) 
nm{K^, A) := {h : h e Hk,, A = s^(h)} . (81) 
Then VA > 0, since 7Yo(^o, A) = 0, we can express the set V. of all possible fading states as 

U Hko^ U |7io(i^o,A)U IJ (7i„(Xo,A)U7Y„(Ko,A)) I. (82) 

l<Ko<M l<Ko<M I l<m<KQ ) 
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Using the defined partitions Hrq, we can rewrite (34) as 

iiiin^i(h) ^Ko=i -^hew^o 

(h)C^(h 



subject to E?n=i ^1 



(83) 



If the optimization is feasible, we have 

M 

e^^ukoIRkoM >r- 



(84) 



J!:o=i 



It is easy to show that we can always achieve (37) although the solution for Tq may not be unique. 

In the following, VKq G [1, M] and Vm e [1, Kq], we drop the dependence of f^, f^, i?^, 
and Sm on h for notational brevity. For Kq e [1, M] and Vh e TY^-q, we let 



f* := [fi*,...,fL] , f := [fi,...,fi^J . 



(85) 



The T* given by the theorem satisfies 



M 



■ Ko 

m=l 



R. 



Ko=l 

From the definition of {H^{Kq, A*)}^Lo and {Hm{KQ, X*)}^=^, we have Vf ^ f* that 

-R 



(86) 



,m=l 



.m=l 



M 



/s:o=i 



■ Ko 
,m=l 



M 
Ko=\ 

M 
Ko=l 

M „ Ko 

= V / Vf,i?,rfF(h) 

M A-,, „ /if,, 

+ [Y^^Ri-Rm] dF{h) 

M Ko ^ / Ko 

+ Jl [Jl^iRi- ^oRm - (1 - r;)Rm-i I dF(h). 

By the convexity of J(i?(h)), we can easily show that 

— < -5-' m = 2,...,Ko. 



(87) 



(88) 
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1) Vh e HoiKQ, A*), since A* < si = Ci/Ri, we have from (88) that A* < Ci/IU, Wi e [1, Kq], 
i.e., Ri < Ci/X*. Then the first summand of the last equaUty in (87) will be 

M „ Ko ^ M „ Ko 

J2 / 5^ r,i?, (iF(h) < - J] / ^f,Q(iF(h). (89) 

2) Vh e HmiKo, A*), m = 1, . . . , /To, since Sm < A* < s^+i, we have 

R-m+l < Rm + (C'm+l — Cm)/X*, Rm-1 < Rm + (Cm-1 " Cm)/X*. (90) 

Using convexity, we have < A* < s^, Vi G [1, m — 1], V/c e [m + 1, Kq]. Hence Vi 7^ m, 
i e [1, ii'o], we have Ri < Rm + (Q — Cm)/X*. Then the second summand in (87) will be 



M i^o /. / Ko 

< 

M Ko „ / Ko \ 

EE/ E ^^(c^ - Cm)/A* rfi^(h) 



M Ko . / Ko \ 

EE/ E + (Q - CJ/A*] +f„i?^-i?„ dF{h) 



^ M Ko r. / Ko \ 

- Y*Y.Y. E^^^'-^- ^^(h)- (91) 

3) Vh e Um{Ko, A*), m e [1, Xq] and since A* = we have Rm-i ^ Rm + (Cm-i - Cm)/X*. 
Since Vi G [1, m — 1], VA; G [m + 1, Kq], we have Sj < A* < s^, from which it follows that 
Ri < Rjn + {Ci — Cm)/\*, Vi 7^ m, m — 1 and i G [1, Kq]. Then the third summand in (87) 
will be 

M K„ / Ko \ 

EE/ E ^^^^ - ^oRm - (1 - T*Q)Rm-l dF{h) 

M Ko ^ I Ko 

^ E E / - E r,[i?^ + (a-cj/A*] 

^f^Rra - r^Rm " (1 " T*)i?.„-l ^ dF{h) 
M Ko ^ / Ko \ 

= EE / - E nid-CJ/X* +{l-T*){Rm-Rm-l)] dF{h) 

^ M Ko I. I Ko \ 

= E E / VL'^^Ci- T*QCm - (1 - T*)C^-l dF(h). (92) 

Xo=i ,„=i'^heH^(ifo,A*) Vi=i / 
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Substituting (89), (91) and (92) into (87), we find 

M 



^heH 



m=l 



-R 



1 ^ 

< -Y 



^ Ko=l t-^heHo(i^o,A*) 



+ E / E - ^0 Cm - (1 - r*)Cm-i dF{h) \ 

^^^Jh&Hm{Ko,X*) Vi=i / J 



.i<ro=i 



Kn 



m=l 



M 



Ko=l 



m=l 



(93) 



Therefore, Vf ^ f*, if r satisfies the total average-rate constraint X]xo=i -^heWxo X]m=i "^m-Rmj > 
^, we have 



E 



hew 



Ko=l 



Ko 



m=l 



M 

> E 



m=l 



(94) 



Hence, r* is the optimal solution to (34) and consequently the corresponding r* and r* are the 
optimal solutions to (31). 



E. Proof of Theorem 4 

Following the Lagrange multiplier method, (40) is equivalent to 

niinr(.),x(.) i?h [Ef=i/ifc^(2'^'=('^) - 1)] - Ef=i A.i?h[rfc(h)r,(h)] 
subject to Eh[Tk{h.)rk{h.)] ^ Rk, A; = 1, . . . , 

Ef=i^.(h) = l, Vh. 
Using the definition of fk{rk(ii-)), we can rewrite (95) as 

minr(.),r(.) Ef=i rfe(h)/fc(rfc(h)) 
subject to Eh[Tk{h)rk{h)] ^ Rk, k^l,...,K, 
ELr,(h) = l, Vh. 



(95) 



(96) 
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Due to the feasibility and convexity of the power region V{f), the original problem (39) has a 
solution. Therefore, given any /j, :— [//i, . . . , /x^]^ > 0, there exists a A* := [A*, . . . , A^]^ > 0, as 
well as optimal rate and time allocation policies r( ) and t( ) for (95). It is clear that we can again 
decompose (95) into two sub-problems. Given A*, we first calculate r^(h) and T^(h) by solving 
(41) and then we determine the water-filling level A* by satisfying individual rate constraints. 

Since /fc(rfe(h)) is convex in rfe(h), it is easy to show that it attains its minimum at rfe,inin(h) = 
logA^ - logi^^l . From the KKT conditions [16, Chapter 5], it follows that for Tl{h) > 0, 
we should have r^(h) = rk^mmi^)- Next, we show that a time allocation strategy t(-) 7^ t*(-), 
excluding the arbitrary time sharing when functions /fe(rfc,inin(h)) have multiple minima, is strictly 
suboptimal. 

Let us first consider the two-user case. 

1) Suppose that for a fading state h, we have /i(ri,inin(h)) < /2(?^2,min(h)) and r(h) is a time 
allocation policy different from T*(h); i.e., Ti(h) = 1 — a and T2(h) = a with a > 0. 

Consider the power cost J(h) := (1 - a)/i(ri,min(h)) + a/2(r2,mm(h)), where r2,min(h) > 
by definition. 

a) If r2,min(h) = 0, then we can let r{(h) = 1, r[{h) = ri,inin(h), and r^(h) = r^(h) = 
such that J'(h) = /i(ri,mm(h)) < -'^(h)- It is clear that T[{h.)r[{h.) > ri(h)ri,min(h), 
and since r2,min(h) = 0, T2(h)r2(h) = T2(h)r2,min(h) = 0. Therefore, if instead of T(h) 
we adopt T'(h), we incur lower power cost without violating the average individual 
rate constraints. 

b) If ?"2,min(h) > 0, we should have /2^''(^2,min(h)) = 0. Let us consider the power cost 
J'(x,h) {1- a + x)/i(ri,i„in(h)) + {a - x)f2 (^r2,min(h)), and define g{x) := 
J(h) - J'(x, h). Since ^(0) = and g^^\0) = /2(r2,n.in(h)) - /i(ri,^in(h)) > 0, there 
exists Aa e (0, a) such that g{Aa) > 0. Therefore, we have r[{h) — 1 — a + Aa, 
r;(h) = ri,inin(h), r^(h) = a- Aa and r^(h) = ^^r2,mm(h), such that J'(Aa,h) < 
J(h). Since r{(h)rl(h) > ri(h)ri,^in(h) and r^(h)r^(h) = r2(h)r2,n.in(h), with T'(h) 
and r'(h) we can afford a smaller power without violating the average individual rate 
constraints. 
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2) For a fading state h satisfying /i(ri,niin(h)) > /2(r2,min(h)) and T(h) ^ T*(h), it can be 
similarly shown that we can have J'(h) < J(h) with alternative resource allocation policies. 
Previous considerations show that the optimal resource allocation policies should follow Theorem 
4 for the two-user case. Similar arguments extend readily to the general i^T-user case as well. 

F. Proof of Theorem 5 

i) From the optimal rate the time allocation policies, we can directly verify the following 
fact. For all k, if the kth component of A increases while other components remain fixed, 
^fc(A) increases whereas Ri{\) decreases for i ^ k. More generally, for any subset /C, if 
we increase for all e and hold the remaining Aj fixed, Ri{X) decreases for i e KP , 
where superscript C here denotes set-complement. 

ii) It can be easily verified that when A = 0, f (A) = 0, and as A t oo, f (A) ] oo. This in turn 
implies that 

a) VA(0) > 0, there exists a < A(0) for which y{cx) < f; 

b) VA(0) > 0, there exists (3 > A(0) for which f (/3) > f. 

iii) Upon defining the mapping A: A(Z)^A(Z-|-1), we can readily verify from Lemma 3 and 
i) that: 

a) The vector A* is the unique fixed point of the mapping T. 

b) The mapping A is order preserving; i.e., A < A' =^ ^(•^) < ^(-^O- 

iv) We now establish one more fact. 

a) If A(0) > A(A(0)) and we define X{1) := A'(A(0)), I = 0,1,..., then X{1) is a 
decreasing sequence and X{1) J, A*. 

b) If A(0) < A(A(0)), then A(/) is an increasing sequence and A(/) ] A*. 
Proof: We verify a) and b) as follows. 

a) By preserving the order, we know that A(/) is decreasing. From ii), there exists a 
a < A(0) for which r(a) < f. In addition, by preserving the order, V/, we have 
X{1) > A^(a). But since r(a) < f, from i), we know that A' (a) is an increasing 
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sequence. Hence, is decreasing and bounded; thus, it must converge to the 

unique fixed point A*, 
b) For A(0) < A(A(0)), we can similarly prove the claim. 
Relying on i)-iv), we are ready to prove the theorem. Notice that ii) guarantees the existence of 
points a(0) and /3(0) such that 

I) a(0) < A(0) < /3(0); II) f{cx{0)) < r; III) r(/3(0)) > r. (97) 

Defining a{l) := A'(a(0)) and f3{l) := A'(/3(0)), we know from iv) that a{l) t A* and 
I3{1) i A*. By preserving the order a{l) < X{1) < (3{l), we have A(/) A*. 



G. Proof of Theorem 6 

For all CSI realizations h, let 
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(98) 
(99) 



and define fk := fk,i. If r is feasible, we have A* for which (59) is satisfied. Let us also 



define Hj := {h : Wk{^)}k=i have J minima}, for J e and users {?j}/=i having the 

"best" channels. Then Vr ^ f*, we have 



K 



fe=l 



Mfc 



1=1 



— Rk 



K 



k=i 



EheHi 



.1=1 



K 



ThiPk,l 



Mi 



1=1 

( Mk 

E ( E 1 + E ~ 

,fe^i, k=l 



1=1 



K 
J=2 



K 



'Mi, 



E K^rk,lpk,l + E K - ^jPij,ll 
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.(100) 



1) Vh G Hi, suppose that user i is selected by f* as the user with the best channel. Since 
l^klkrJhk < K < l^klk,ii+i/hk, k e [l,K], we have ^"''^ ^"''^"^ < At < ^"''^ 



Pk,ii-Pk,ii-i 

k k 



Pk,l*+1-Pk,l* 



Therefore, 

< Pk,i*k + (^Mfc-i ~ Ck,i{)/K^ pk,ii+i < Pk,i* + {Ck,i*+i - Ck,i*)/x*k. (101) 
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By the convexity of Tfe(x), we have Hklk,i/hk < K < IJ'klk,i' / hk, V/ e [1,/^ - 1], V/' e 
[II + 1, Mfe]. Hence, VZ II, I e [1, M^], we have p^,, < pk,i* + {Ck,i - Ck,i*)/Xl. Then the 
first sum of the last equality in (100) will be 
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(102) 
(103) 



— -E'hGWi ^ ^ T'k,lCk,l — Ci^i 
_k=l 1=1 

where we used the allocation policies V/c ^ i, Ck,i* — XkPk,ii > ~ KPhi* obtain 
inequality (102). 

2) Vh e Hj, J > 1, suppose that users ij are selected by f* as the users with best channels, and 
define C^in := Q.^i*^ - KjPij,i*.' J = '^,---,J- Noticing that V/c ^ ij, Ck,i* - X*kPk,ii > C^n 
and that VZ II, I e [1, Mk], Pk,i < Pk,ii + {Ck,i - Ck,ii)/\l, we have 
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K Mk 
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k=l 1=1 j=l 
K Mk J 

,k=l 1=1 3=1 
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_k=i 1=1 

Substituting (103) and (104) into (100), we have 
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(105) 



Therefore, Vr ^ f*, if r satisfies individual rate constraints £'h Tfe,/pfc,;j > Rk, we have 



■ K Mk 

.fe=i 1=1 



> E. 



■ K Mk 

E E^M^fc-' 
.fe=i j=i 



(106) 



Hence, r* is the optimal solution to (56) and consequently the corresponding r* and r* are the 
optimal solutions to (55). 

Similar to Lemma 3, we can show that A* is almost surely unique. Define Q as the set of all 
feasible rate vectors. If f is feasible, there must be a boundary point Vg of Q for which f < fg. Let 
\g denote the Lagrange multiplier corresponding to Vg. Then for A(0) < \g, there exists /3 > A(0) 
for which r(/3) > r. With this guarantee replacing the counterpart ii)-b) and following the lines 
in the proof of Theorem 5 (Appendix F), we can show that A* can be iteratively computed by 
Algorithm 4 for any positive initialization A(0) < A^. 
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rate reward 



Fig. 1. The functions fi{R), f2[R) and J{R) when wi < TO2 and -^^^ < -^^^ (the dash-dotted curve for J{R) is slightly 
offset for easy visualization). 




rate reward 



Fig. 2. Convex envelope J{R) in the finite-AMC-mode case (the dash-dotted curve for J{R) is slightly offset for easy 
visualization). 
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3. Power regions for the infinite-codebook case when two users have identical SNRs: h\/{NoB) = h2/{NoB) — dBW. 
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4. Power regions for the finite-AMC-mode case when two users have identical SNRs: hi/{NoB) — h2/{NQB) = dBW. 
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Fig. 5. Power regions for the infinite-codebook case whien two users have 10 dB difference in SNRs: hi/{NoB) = 10 dBW, 
and h2/{NoB) = dBW. 
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Fig. 6. Power regions for the finite-AMC-mode case when two users have 10 dB difference in SNRs: h\/[NoB) = 10 dBW, 
and h2/{NoB) = dBW. 
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-A- optimal policy under weighted sum rate constraint over policy A 
-B- optimal policy under weighted sum rate constraint over policy B 
optimal policy under individual rate constraints over policy A 
optimal policy under individual rate constraints over policy B 




Fig. 7. Power savings for the infinite-codebook case when two users have identical SNRs: hi/{NoB) = h2/{NoB) = dBW. 
(Pohcy A: equal time allocation and separate water-filling; Policy B: equal time allocation among users and equal power per fading 
state for each user.) 
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-B- optimal policy under weighted sum rate constraint over policy B 
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Fig. 8. Power savings for the finite-AMC-mode case when two users have identical SNRs: hi/{NoB) = h2/{NoB) = dBW. 
(Policy A: equal time allocation and separate water-filling; Policy B: equal time allocation among users and equal power per fading 
state for each user.) 
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Fig. 9. Power savings for the infinite-codebook case when two users have 10 dB difference in SNRs: hi/{NoB) — 10 dBW, 
and h2/{NoB) — dBW. (Policy A: equal time allocation and separate water-filling; Policy B: equal time allocation among users 
and equal power per fading state for each user.) 



45 



